Extended Hubbard model: Charge Ordering and Wigner-Mott transition 
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Strong correlations effects, which are often associated to the approach to a Mott insulating state, 
in some cases may be observed even far from half-filling. This typically happens whenever the 
inter-site Coulomb repulsion induces a tendency towards charge ordering, an effect that confines the 
electrons, and in turn favors local moment formation, i.e. Mott localization. A distinct intermediate 
regime then emerges as a precursor of such a Wigner-Mott transition, which is characterized by both 
charge and spin correlations, displaying large mass enhancements and strong renormalizations of 
other Fermi liquid parameters. Here we present a careful study of a quarter filled extended Hubbard 
model - a simple example where such physics can be studied in detail, and discuss its relevance for 
the understanding of the phenomenology of low-density two dimensional electron gases. 



I. INTRODUCTION 

Early theoretical and experimental investigations of 
two dimensional electron gases (2DEG) largely focused 
on disorder effects that dominate the so-called "diffu- 
sive" regime ksT -C H/t, t being the impurity scat- 
tering rate. In this situation, which is best realized in 
relatively low-mobility materials, the interaction effects 
are expected to be weak and coherent multiple-scattering 
processes dominate - ultimately leading to the forma- 
tion of bound electronic states through Anderson local- 
ization. This impurity effect was first predicted in the 
famed scaling theory of localization, 1 for non- interacting 
and weakly disordered 2DEG systems, and was quickly 
extended to include the weak interaction corrections. 2 
The predicted logarithmic rise of the resistivity at low 
temperature was soon confirmed by experiments on thin 
metallic films and two-dimensional semiconducting sur- 
faces. 3,4 All this intensive activity around the 2DEG con- 
tributed to the emergence of a widely held opinion that, 
in these systems, even a minute amount of (impurity) 
disorder will localize all the electronic states at T = 0. 
If this were true, then there should not exist a sharp 
metal-insualtor transition (MIT) in any two dimensional 
system; the density and/or temperature dependence of 
transport should simply reveal a gradual crossover from 
weak to strong localization. 

Given this conventional lore focusing on disorder ef- 
fects, the 1994 pioneering experiment of Kravchenko et 
al. 5 provided quite a surprise. It reported low tempera- 
ture transport behavior in ultra clean 2D Silicon samples 
which, around this time, became available due to tech- 
nological advances in fabricating semiconductor devices. 
Kravchenko's work presented first evidence for dramatic 
changes in the temperature dependence of the resistivity 
in a narrow density range, suggesting the possibility of a 



metal-insulator transition (MIT) in the 2DEG. This re- 
sult - quite surprisingly - passed almost unnoticed for sev- 
eral years, until confirmed by independent transport and 
magnetic response measurements on clean Si-MOSFET 6 , 
as well as on other semiconductor heterostructures (e.g. 
GaAs/AlGaAs) 7 . These results have attracted a great 
deal of attention to this field, triggering a surprisig va- 
riety of proposed theoretical scenarios 8 , many of which 
were quickly ruled out on experimental grounds 9 . 



A. Strong correlations revealed 

The next ten years produced significant new infor- 
mation from complementary studies by several groups, 
which reported singular enhancement of the effective 
mass, 9 to* ~ (n — n c )~ l , and emphasized the central 
role played by spin physics in the low-density 2DEG sys- 
tems. In particular, experimental measurements of the 
electronic spin susceptibility in the proximity of the MIT 
shown a Curie- Weiss behavior xl n — 91^%/T, suggesting 
almost total conversion of the electrons into local mag- 
netic moments below the critical density. 10 ~ 12 While the 
inital attention concentrated on the role of disorder 1,3 , 
more recent experiments 9 made it increasingly clear that 
many key experimental features seem to persist even in 
the cleanest samples, thus pointing to effects intrinsic 
to the low-density 2DEG. All these experimental results 
have contributed to validate the emerging idea 14,15 that 
strong interactions are the primary driving force behind 
the instability of the Fermi liquid phase in favor of an 
insulating state in the low density 2DEG. 

In very recent work, 1 ' 15 radically new ideas have been 
put forward to explain the observed behavior, and to rec- 
oncile the early viewpoint of Wigner with the possibility 
of Mott physics in the low density regime. According to 



this Wigner-Mott scenario 14 ' 15 the rich phenomenology 
of the two-dimensional electron gas (2DEG) is ultimately 
determined by the competition between the long-range 
Coulomb interaction and the kinetic energy. As first 
noted by Wigner, 16 ' 17 both the the Coulomb interaction 
Ec and the kinetic (Fermi) energy Ep decrease with de- 
creasing density n, but their ratio r s = Ec/Ep, i.e. the 
Wigncr-Seitz radius, increases. Thus, the Coulomb re- 
pulsion is expected to dominate at very low densities 
(r s 3> 1) and low enough temperatures, leading to the 
formation of a charge-ordered state - a Wigner crystal. 
In this regime each electron forms a bound state within 
a potential well formed by the repulsion from other elec- 
trons, and must overcome an activation gap to escape. In 
the opposite regime (r s <C 1) the quantum fluctuations 
lead to the eventual melting of the Wigner lattice, and 
the consequent formation of a homogeneous Fermi liquid 
state. In between these two limits the system is expected 
to undergo a series of non-trivial transformations, mak- 
ing the description of the phase diagram of the 2DEG 
one of the most challenging tasks of modern condensed 
matter physics. 



B. Driving force: interactions or disorder? 

This Wigner-Mott scenario should be contrasted with 
the alternative physical picture as first proposed in the 
early theoretical works of Refs. [18,19]. Here disorder is 
envisioned as the primary driving force for electron local- 
ization, and the insulator consists of Anderson-localized 
electrons bound to impurities. In this scenario, the pri- 
mary role of the interactions is to stabilize the metallic 
state at higher densities, a possible mechanism that has 
been explored in relatively recent theoretical work. 13 The 
key theoretical challenge, therefore, is to carefully iden- 
tify the precise consequences of each of these physical 
pictures, and assess their respective relevance in light of 
experiments. In our work we deliberately ignore all dis- 
order effects, and focus on investigating the predictions 
of the Wigner-Mott picture, which we believe can acount 
for most qualitative aspects of the puzzling experimental 
features. 

To understand how the Mott physics, usually ascribed 
to half-filled narrow bands, emerges out of the low- 
density 2D systems, one should consider the following ar- 
gument. A barely melted Wigner crystal is dominated by 
the local spin correlations while the short-range Coulomb 
repulsion largely precludes double-occupations. There- 
fore, it becomes reasonable to view a metal in the vicin- 
ity of Wigner crystalization as a system on the brink of 
Mott localization. Physically, in the r s ^> 1 regime, the 
Coulomb interaction proves so strong that it keeps all 
the electrons at bay, even after the crystal has melted. 
The volume fraction available for each electron is thus 
significantly reduced, giving rise to the "confinement" of 
each electron. This situation is locally close to half-filling 
and the Mott regime, hence providing an explanation for 



the (observed) similarities of a dilute 2D electrons gas 
with conventional Mott materials characterized by nar- 
row half-filled bands. 

While of plausible significance for diluted 2DEG, 
the phenomenon of Wigner-Mott localization applies to 
many more systems. In fact, based on the general argu- 
ments of Wigner, one may expect that it will emerge for 
any partially filled band, where inter site (long-range) 
Coulomb interactions are able to induce charge order- 
ing. These effects have generally been little explored, 
and were so far mostly studied in (quasi) one-dimensional 
systems. 20,21 However, recent detailed experimental stud- 
ies on charge ordering phenomena in layered organic 
molecular crystals, 22 and on the cobalt oxide Naa;Co02, 23 
have triggered several theoretical investigations of ex- 
tended Hubbard models which includes both on-site and 
nearest neighbor interaction. 24 ~ 26 Some aspects of the re- 
sulting behavior may be influenced by material-specific 
details, such as the form of the lattice. In compounds 
with triangular lattice structure, for example, geometri- 
cal frustration leads to the competition of several ground 
states 27 . Although still relatively few, all these exam- 
ples indicate that Wigner-Mott localization is clearly a 
very general phenomenon, which is yet to be explored in 
detail. Its experimental realizations are often found in 
two-dimensional systems, where device geometries (e.g. 
gating the 2DEG) allows easy control of carrier density, 
which facilitates accessing the low density transition re- 
gion. Still, as a matter of principle, this phenomenon is 
not restricted to d = 2, and thus should be found in all 
dimensions. In addition, the charge density wave order- 
ing (CDW) underpinning can emerge for both long-range 
inter-site (Coulomb) interactions, but also for sufficiently 
strong short range (e.g. nearest-neighbor) repulsion. 

To gain insight into its generic features, in this paper 
we restrict our attention to the simplest model that il- 
lustrates its fundamental mechanism: Mott localization 
as driven by charge ordering . This is accomplished by 
carefully examining the quarter-filled extended Hubbard 
model (EHM), 28 ' 29 which we solve using single-site dy- 
namical mean field theory 30 (DMFT) methods and a 
combination of several quantum impurity solvers. The 
EHM has already been studied in some detail in previous 
work by Pietig et al. 31 , using the DMFT approximation. 
These authors found a charge-ordered phase that forms 
at large values of the inter-site interaction, by using a 
combination of non-crossing approximation (NCA), ex- 
act diagonalization (ED), and numerical renormalization 
group (NRG) calculations. This work demonstrated the 
increased importance of correlation effects in the charge- 
density wave (CDW) phase, leading to the formation of 
a strongly renormalized quasiparticlc with the density of 
states displaying a narrow peak at the Fermi level. How- 
ever, the complete V-U phase-diagram for this model 
has not been determined and, more importantly, no evi- 
dence for the formation of a Wigner-Mott insulating state 
has been reported. These and some other interesting as- 
pects of the 2DEG have been clarified in a more recent 



short paper, 15 where the DMFT equations for the EHM 
have been solved at low but finite temperature using 
the recently-developed continuous-time Quantum Monte 
Carlo algorithm. 32,33 This work identified the existence 
of a strongly correlated metallic CDW phase, separating 
the homogeneous metal from a Wigner-Mott insulating 
phase present for larger values of the non-local interaction 
and small enough values of the short-range interaction. 
In contrast, a direct transition between the Fermi liquid 
metal and the Wigner-Mott insulating state was found 
at larger values of the local correlation, triggered by the 
tendency to charge ordering. 

Despite this progress, however, several important phys- 
ical issues remain to be clarified, in order to correctly 
understand the phenomenology of the Wigner-Mott sce- 
nario and its possible relevance for the 2DEG problem. In 
particular, even within the considered DMFT solution of 
the EHM, the following questions need to be addressed: 

1. What are the conditions required to find an inter- 
mediate strongly correlated CDW metallic (CDW- 
M) phase? 

2. What is the temperature dependence of this phase? 

3. What is the nature of the phase transitions between 
the different phases and how does the critical be- 
havior depend on temperature? 

4. Is there a regime of phase coexistence between the 
metal and the insulating state, allowing the possi- 
bility of phase separation? 

This paper presents careful calculations providing precise 
and convincing answers to all these questions. In partic- 
ular, we show the existence of an intermediate CDW-M 
phase, that separate the homogeneous metal from the 
Wigner-Mott insulator, for all temperatures smaller than 
a critical value T = T c . The shrinking of the CDW- 
M phase corresponds to temperature dependent behav- 
ior of the Wigner-Mott transition, similar to the Pomcr- 
anchuk effect. 34 We identify two transitions in the phase- 
diagram of the model, namely a charge-ordering transi- 
tion from the homogeneous metal to CDW-metallic state 
at V = V c i(T) and a Wigner-Mott MIT taking place at 
V = V C 2(T). We demonstrate the continuous character 
of the MIT for T < T c and we give clear indications 
for the existence of phase-coexistence, i.e. a sharp first- 
order transition, at higher temperatures. Furthermore, 
we present a detailed study of the evolution of the two 
critical lines V c \ and V c i as a function of increasing tem- 
perature from T = 0. 

The rest of the paper is organized as follows. In section 
II we introduce the main model and the relative notation, 
together with a brief overview of the (numerical) methods 
used in this work. In section III we present the results 
concerning the Wigner-Mott transition in the EHM and 
we discuss the related phase-diagram. We conclude in 
section IV, where we present our perspective on the rele- 
vance of our picture for the 2DEG and discuss its relation 
to alternative theoretical scenarios. 



II. MODEL AND METHODS 
A. Extended Hubbard model 

The extended Hubbard model is the simplest model 
that capture the interplay between strong correlation and 
charge-ordering effects. This model includes both a local 
Coulomb interaction, represented by the familiar Hub- 
bard U term, and a non-local (inter-site) repulsion V . 
While the presence of this inter-site repulsion V may 
induce the formation of a charge-ordered phase, it also 
proves capable to enhance the effectiveness of the on- 
site repulsion U, leading to the formation of a strongly 
correlated physics even away from integer-filling. The 
lattice structure of the model is introduced in order to 
capture the crystalline order of the 2DEG. The lattice 
constant is then constrained by the requirement that each 
cell contains two lattice sites, which can be regarded as 
precursors of the interstitial and vacancies in the Wigncr 
crystal phase, corresponds to an area of nr 2 a B . The cor- 
responding extended Hubbard model Hamiltonian takes 
the form: 
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Here, U and V are the on-site and the nearest neighbor 
interactions, respectively. The parameter t is the hop- 
ping amplitude, c ia {ci a ) arc the creation (annihilation) 
operators and n, = n,-f + n^ is the occupation number 
operator on site i. The chemical potential fi is adjusted to 
enforce the quarter filling constraint. To this end we can 
restrict ourselves to consider a bipartite lattice structure 
with sublattices A and B. 

We solve the problem posed by the Hamiltonian Eq. 1 
using single-site DMFT method, 30 which amounts to 
mapping the previous lattice problem onto that of a sin- 
gle impurity, coupled to an effective bath, to be self- 
consistently determined. This approximation becomes 
exact in the infinite coordination number limit, 35,36 pro- 
vided the hopping parameter and the inter-site interac- 
tion are rescaled as t — > t/y/z and V — > 2V/z. In this 
limit, the interaction between electrons on neighboring 
sites i and j are treated in the Hartree approximation, 
and the Hamiltonian Eq. 1 assumes the form: 

H = -t ^2 ( c ^^ c h,a + h.c) -iJ,y^ni + Uy~] ntfriii 
+ v Yl i( n i) n o + n i(nj) - (ni){nj)) 

i£A,j£B 

For simplicity, the terms linear in the occupation opera- 
tors can be easily absorbed in a redefinition of the chem- 
ical potential /ia = M — 2V{ub)^ /J-b = I 1 ~ %V(tia}- The 
Weiss field, describing the properties of the effective bath 
in the DMFT approximation, reads: 

Go-A.B^n) = iu n - HA,B - &A,B(iUn), 



(2) 



where A^ j 5(iw„) is the hybridization function of the as- 
sociated effective single impurity Anderson model. The 
self-consistency condition for the model at hand assume 
the following expression 30 : 



Ga,b = (b,a / ds 
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with Ca,b = iu n + p A ,B - Z AjB (iu) n ), Y> a ,b being the 
local self-energy function and po(e) = X)k^( £ — e s )) 
the non-interacting density of states corresponding to 
the chosen lattice. In the following we focus on the 
simple semi-circular model density of states, for which 
po(e) = VD 2 - e 2 /4nt 2 with D = 2t = 1 fixing the en- 
ergy units of the problem. The self-consistency equations 
then reduce to: 
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B. Methods 

In this section we briefly review the methods of solu- 
tions used in this work, emphasizing the advantages and 
the drawbacks for each of them. More detailed review 
of the same methods can be found elsewhere in the lit- 
erature, cf. Ref. 30,32,37. We begin our discussion with 
the simplest used impurity solver, namely the four slave- 
boson (SB4) method of Kotliar and Ruckenstein. 38 At 
the mean field level (saddle-point), this method at T = 
proves to be equivalent to the well-known Gutzwiller vari- 
ational approximation, but it also allows for extensions 
to T > 0. Within the SB4 method the effective impurity 
problem is reduced to the solution of a set of algebraic 
equations for the appropriate slave-boson variational pa- 
rameters. Solving these equations is a relatively simple 
task, which allows for quite an accurate numerical so- 
lution, and a very precise characterization of the lead- 
ing critical behavior. Although the SB4 method varia- 
tionally calculates the Fermi liquid parameters specifying 
the quasiparticles and by neglecting incoherent parts of 
the spectrum, this method is known to generally provide 
qualitatively and largely even quantitatively accurate so- 
lutions when applied within the framework of DMFT 
theories. In the SB4 approach, the low energy part of 
the local Green's function is parametrized in terms of 
the quasiparticlc weight Zi at site i. In the CDW phase 
these quantities differ on the two sublattices i = A, B. 
The quasiparticlc weight Zi[ei,di\ is expressed through 
the mean field slave-boson parameters ef and d 2 , which 
are equal to the probability that the site is empty and 
doubly occupied, respectively. These parameters and ul- 
timately the impurity Green's functions GfAjB, for ei- 
ther the A and the B sublattices, are determined by con- 
strained minimization of a suitable local free energy Ti 



parametrized by a Lagrange's multiplier Xf. 
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with p.A,B = I 1 ~ 2V n B.A — ^a.b- The (low energy part 
of the) local Green's functions are determined by the 
relation Ga,b = Za,bGjajb- The DMFT equations 
are closed by supplementig these relation with the self- 
consistency condition 3. The chemical potential p and 
occupation numbers ua and ub also have to be calculated 
self-consistcntly in order to enforce quarter-filling con- 
straint n A + n B = 1, where n A ,B = f XL„ GfA,B(iw n )- 

To supplement the SB4 results with more accurate 
methods we solved the DMFT equation for the EHM 
with (numerically) exact methods. In particular in the 
T = limit we implemented a Density Matrix Rcnor- 
malization Group (DMRG) impurity solver. The DMRG 
idea have first been introduced in 1992 by S. White 39 to 
deal with one-dimensional quantum lattice problems, 40 
for which the Wilson's real-space blocking scheme was 
observed to fail in describing the correct solution. In 
more recent work, 37,41 DMRG has been adapted to solve 
the self-consistent DMFT equations. The algorithm is 
based on a recursive exact diagonalization (ED) of the 
associated quantum impurity problem with an increasing 
size of the effective bath. Similarly to other renormaliza- 
tion group based method, the DMRG provides a suitable 
method to restrict the solution of the effective problem 
to its relevant subspace. This is achieved using a clever 
method for limiting the exponential growth of the Hilbert 
space dimension, based on the analysis of the reduced 
density matrix. Starting with a small effective bath in 
the form of two linear chains, 30 the effective impurity 
problem is exactly solved, and then the size of the bath 
is recursively increased. At each step the Hilbert space of 
the problem is constructed using the basis formed by the 
first M (i.e. the "most probable") eigenvectors of the re- 
duced density matrix p = Tr| B?iu \gs)(gs\, where \gs) is the 
(approximated) ground state and the environment is rep- 
resented by the impurity plus one of the two bath' chains. 
For every fixed size of the impurity problem, the impu- 
rity Green's function are evaluated using successive ap- 
plications of the Lanczos' method. 42 The self-consistency 
condition 3 is then used to update the Hamiltonian pa- 
rameters of the problem. 

The recursive nature of the DMRG method permits a 
systematic improvement of the quality of the solution. 
In addition, this method has the advantage of treating 
on the same footing both low- and high-energy scales, in 
contrast to standard RG technique, which principally fo- 
cus on the low-energy spectra of the problem. However, 
the particular one-dimensional topology imposed by the 
DMRG algorithm makes this method suffering of the fi- 
nite size effects. In practice effective baths as large as 
L = 20 sites have to be reached in order to obtain a 
satisfying solution of the model. 



While the DMRG algorithm it is straightforward at 
zero temperature, its extensions to finite temperature 
regime are non-trivial. To this end, we have solved the 
DMFT equations using alternative algorithms. The first 
finite temperature method is (complete) Exact Diago- 
nalization of the discretized effective impurity problem. 
This method permits to access the full spectrum of the 
problem, thus allowing for the calculation of the exact 
impurity Green's function and other observablcs. The 
Hamiltonian parameters are self-consistently determined 
using an adaptive method that permits to minimize the 
(in principle large) finite size effects. 43 This method is in- 
deed known to have an almost constant scaling with the 
bath' size. On the other hand, a small number of degrees 
of freedom in the problem prevents the access to the very 
low energy physics of the model, thus a significant num- 
ber of bath' sites (N s > 5) becomes necessary to have a 
good description of the model solution. Our calculation 
has been performed using N s = 7, which is the largest 
accessible number of sites for this method. 

The second method we used to solve the DMFT equa- 
tions at finite temperature is the Continuous Time Quan- 
tum Monte Carlo, in the implementation of Ref. [32]. 
The CTQMC is a statistically exact method based on the 
Monte Carlo samplings of the diagrams obtained by per- 
turbativc expansion of the impurity problem with respect 
to the hybridization. The CTQMC algorithm of Ref. [32] 
has proven to be a highly reliable and stable method, 
that permit to capture the low temperature regime of 
the model. On the other hand, the statistical nature of 
the method put some, indeed non-stringent, limitations 
to its applicability in the proximity of a phase transition. 
In this regime, the solution requires a very large number 
of Monte Carlo samplings and a large number of itera- 
tions of the self-consistency algorithm to overcome the 
enhanced fluctuations and the critical slowing down of 
the solution. 
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Figure 1: (Color online) Non-interacting DOS of sublattice A 
(solid line) and B (dashed line) for 5 ~ 0.2. The system is in 
a CDW-M state and correspondingly the DOS is split in two 
lobes. The Fermi energy, corresponding to ui = 0, falls inside 
the lower lobe as results of the quarter filling constraint. 



occupation imbalance between the two sublattices. The 
density of states is non-zero only in the energy interval 
defined by the 1 — j3/a < condition. The system un- 
dergoes a charge ordering at V = V c with the formation 
of a spectral gap A = VS (for V > V c ). In this regime 
the original quarter-filled band splits into two bands, but 
the Fermi level remains inside the lowest band, which re- 
mains half- filled, cf. Fig. 1. Thus a spectral gap at the 
Fermi level can only be opened by increasing the local 
correlation U. 

In this limit we can think at the system as consisting of 
one sublattice A with occupation close to half filling, and 
the other B as being nearly empty. An estimate of the 
effective hopping between sites in the sublattice A can be 
obtained integrating out (virtual) high energy processes, 
corresponding to hops through the sites of sublattice B. 
At V ^> i, to leading order the effective hopping between 
the sites of the A sublattice then takes the form t e ff = 
t 2 /V. 



In the non-interacting limit U — the problem posed 
by the Hamiltonian Eq. 2 can be easily solved analyti- 
cally. In particular here we are interested in obtaining 
the analytic solution corresponding to the DMFT ap- 
proximation, in order to get some physical insight from 
the model. In the non-interacting limit the DMFT equa- 
tions can be recast in the form: 



t 2 Gi 



G/iz) =z + ha -t 2 G B = 

Gg 1 (z) = z + /i B - t 2 G A = /3- t 2 G A 

and solved for, say, Ga giving: 
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All the solutions to this equation can be expressed in 
terms of a single parameter <5 = ha — "s representing the 



III. THE WIGNER-MOTT TRANSITION 
A. Phase-diagram 

In this section we present results concerning the phase- 
diagram of the EHM as a function of the local correlation 
U and the inter-site interaction V. In the following we 
use different methods to investigate the various phases of 
the system as a function of the temperature T. 

Our results show the existence of three qualitatively 
different phases, which can be classified by two order pa- 
rameters: i) the CDW staggered density 6 = ua — ub and 
ii) the quasi-particle residue (i.e. renormalization con- 
stant) Za,b at each sublattice. The different phases of 
the model are characterized as following: 

• Fermi liquid. A featureless homogeneous metallic 
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Figure 2: (Color online) DMFT phase diagram of the ex- 
tended Hubbard model obtained with a slave-boson (SB4) 
impurity solver. The figure shows the phase-boundaries ob- 
tained for different values of the temperature T. At T = 
(solid line) the solution shows the existence of a large CDW 
metallic phase at all values of the local correlation U. This 
region narrows as the temperature T is increased from zero 
to T = 0.01 (dashed line) to T = 0.05 (dot-dashed line), but 
it never disappears. 



state, corresponding to ua = ns = 1/2 and Za = 
Z B >0. 

CDW-M. A CDW-ordered metallic state with 
Ua > ns and Zb > Za- The sublatticc A has 
an occupation near to the half-filling condition 
(riA — 1), thus closer to Mott localization. As a 
consequence of this the quasi-particles weight Za is 
substantially reduced, corresponding to a enhance- 



ment of their effect mass 



JrriA 



• CDW-I. A Wigner-Mott insulating phase. In this 
regime ha 3> ns but while Zb ~ 1, sublattice A 
has an identically zero quasi-particle weight Za = 
0. This corresponds to the Mott localization of the 
electrons on the "nearly half-filled" sublattice A. 

The U — V phase diagram as obtained from our SB4 
solution for increasing values of the temperature T is pre- 
sented in Fig. 2. At zero temperature, we obtain a contin- 
uous transition from the homogeneous metallic phase to 
the CDW- metal. This is followed by a second continuous 
transition to the CDW-insulator as the inter-site inter- 
action V is further increased. The intermediate CDW 
metallic phase, while decreasing in size at large U, re- 
mains of finite extent even in the limit U — > oo limit. 
Upon increasing the temperature T, we observe a shrink- 
ing of the intermediate CDW-M phase at large values of 
U. Still, no signs for the disappearance of this phase have 
been observed with the slave-bosons methods. This re- 
sult is in contrast to that initially obtained in Ref. 15 at 
very low but finite temperature T = 0.01 using a more 
sophisticated method. Despite the qualitative validity of 
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Figure 3: (Color online) DMFT phase diagram of the ex- 
tended Hubbard model obtained with CTQMC method. The 
results illustrate the evolution of the phase-boundaries for the 
Mott-Wigner transition as a function of increasing tempera- 
ture. For the lowest value of the temperature T — 0.0025 
(solid line) we observe the well established presence of a 
CDW-M separating the homogeneous metal from the CDW- 
insulator. This phase-region narrows upon increasing the 
value of the temperature to T — 0.0033 (dashed line) and 
T = 0.005 (dotted line), till a closure of the phase-boundaries 
is observed for T ~ 0.01 (dot-dashed line) at U ~ 2. 



the SB4 method, this result clearly puts in question the 
main features of the phase-diagram for the EHM. In par- 
ticular, it suggests that the intermediate CDW-M phase 
at intermediate to large U is very fragile to tempera- 
ture. This observation calls for a more detailed study of 
the temperature dependence of the metal-insulator phase 
boundary, which we examine in the following sections. 

The results from CTQMC solution of the DMFT equa- 
tions for the EHM arc summarized in the phase-diagram 
in Fig. 3. While this method cannot capture the zero 
temperature properties of the model, it provides an excel- 
lent and reliable solution at very-low temperature scales. 
At the lowest accessible temperature T = 0.0025 we ob- 
serve a continuous transition between a homogeneous 
metallic state to a CDW-metal for low values of V and 
every value of U . At this value of the temperature the 
CDW-M phase is found to be destabilized towards a 
CDW-insulating phase upon increasing the inter-site in- 
teraction V . Upon increasing the temperature we ob- 
served a narrowing of the intermediate CDW-M in favor 
of the Wigner-Mott insulating state, until a merge of the 
boundary lines is observed for temperatures of the order 
T ~ 0.01 (that is two orders of magnitude smaller than 
the bare bandwidth). 

To better clarify the issue of closing boundary-lines 
at low temperature we have solved the DMFT equa- 
tions with an exact method at zero temperature, i.e. the 
DMRG. Furthermore, we have complemented this inves- 
tigation with a powerful ED solver, in order to access the 
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Figure 4: (Color online) DMFT phase diagram of the ex- 
tended Hubbard model obtained with DMRG and Exact Di- 
agonalization (finite temperature). The figure shows the evo- 
lution of the phase-boundaries as a function of increasing tem- 
perature from T = (solid line) to T — 0.0033 (dashed line), 
to T = 0.01 (dot-dashed line). The intermediate CDW-M 
observed at T = is found to narrows upon increasing the 
temperature, until a closure of the boundary-lines is obtained 
around U ~ 2 and for T ~ 0.01. The DMRG results have been 
obtained with a bath of L = 30 sites, whereas the ED calcu- 
lation have been performed with a number of sites N — 7. 



finite temperature properties of the system. The ED- 
DMRG results for the extended Hubbard model problem 
are condensed in the U — V phase diagram presented in 
Fig. 4. The phase-diagram is in very good agreement 
with that obtained from the CTQMC solution of the 
model, already presented in Fig. 3. Our calculations con- 
firm the existence of an homogeneous metallic phase at 
low values of V and any value of the local correlation 
U. The stability of the homogeneous metallic phase is 
related to the quarter filling condition, making the local 
interaction nearly ineffective at small V. For any fixed 
value of U and for V larger than a critical value V c \ 
almost independent of the temperature T, the system 
shows a continuous phase-transition towards a CDW-M 
state. The evolution of the two sublattices upon fur- 
ther increasing V is different. One sublattice gets nearly 
empty and becomes, for sufficiently large U, a band in- 
sulator. The other sublattice, with a filling closer to 1, 
undergoes a continuous Mott transition towards a CDW- 
Insulating state for U large enough and V larger than a 
critical value V C 2 weakly depending on the temperature 
T. In particular the zero temperature solution of the 
model shows the persistence of the intermediate CDW- 
M for any value of the local correlation U. This phase 
is observed to narrow upon increasing the temperature 
until a closure of the boundary lines at U ~ 2 is obtained 
for T ~ 0.01, in excellent agreement with the CTQMC 
solution of the model. 
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Figure 5: (Color online) Order parameter S — ua — fiB as a 
function of the inter-site interaction V, for U — 2 and increas- 
ing values of the temperature T. The metallic CDW phase 
appears for V > V c i — 0.84, while a CDW-insulator forms for 
V > V c2 ~ 1.35. Within SB4 solution, the T - transition 
is continuous, while for T > a small discontinuity in CDW 
order parameter is observed, that is a possible artifact of the 
method. 



B. Generalized Pomeranchuk effect 

Next, we would like to discuss the physical implica- 
tions of the observed temperature evolution of the phase 
diagram, in relation to the interesting non-monotonic re- 
sistivity behavior observed in very clean 2DEG samples. 3 
As a general trend we observed a shrinking of the CDW- 
M region as temperature is increased from T = 0. This 
effect is most pronounced at the CDW-M to CDW-I 
boundary V c 2 (U), i.e. in the region where a heavy Fermi 
liquid forms. From the physical point of view, this tem- 
perature dependent behavior reflects the entropy gained 
in destroying the (heavy) Fermi liquid to form localized 
magnetic moments in the Wigner-Mott (CDW-I) phase. 
This mechanism of entropy release is indeed similar to 
that anticipated in the early work of Pomeranchuk, 44 who 
speculated about the general problem finite-temperature 
solidification of 3 He. 34 

It should be emphasized, however, that this "entropic" 
destruction of a strongly correlated Fermi liquid is a very 
general phenomenon, which does not necessarily require 
a first order phase transition, as postulated by Pomer- 
anchuk, or a phase coexistence with separation of the 
phases. 45 ~ 48 For example, the same effect can be observed 
in many heavy fermion systems, where the destruction of 
the heavy Fermi liquid for temperatures larger than the 
"coherence" temperature T* ~ Tjc, coincides with a re- 
lease of extensive spin entropy S ~ fcs In 2, corresponding 
to the localization of the /-electrons. The corrrspond- 
ing resistivity maximum represents a temperature-driven 
crossover rather than a sharp a phase-transition, related 
to the fact that the two phases do not have to differ by 
symmetry. In the Wigner-Mott picture of the 2DEG, we 
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Figure 6: (Color online) The two panels show the behavior of 
the order parameter S = ua ~ na as a function of inter-site 
interaction V, for (7 = 2 and different values of the tem- 
perature T. (a) Results from zero temperature calculations 
with DMRG method, (b) CTQMC solution for several tem- 
peratures T = 0.01 (black dots), 0.005 (red squares), 0.0033 
(green diamond), 0.0025 (blue triangles). Both T = and 
finite temperature calculations display a rather conventional 
(V — Vd) 1 ' 2 critical behavior near the CDW transition, sim- 
ilar to that observed in the non-interacting limit. 



thus recognize in the entropy of the high temperature 
state as the driving force for electron localization. 



C. Charge ordering 

To better understand the nature of the continuous 
transitions observed in the phase-diagrams of Sec. II, we 
now investigate the behavior of the charge order parame- 
ters 8 = n^ — riB as a function of V for different values of 
the temperature. The symmetry breaking associated to 
this transition corresponds to the occupation unbalance 
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Figure 7: (Color online) Quasi-particle weight Za.b as a func- 
tion of V for different temperatures and U — 2. At T = 
(solid line) we observe a linear vanishing of the non-empty 
sublattice renormalization constant, while the constant cor- 
responding to nearly-empty sublattice flows to one. At larger 
value of the temperature T = 0.01 (dot-dashed line) and 
T — 0.05 (dashed line) the SB4 predicts a small jump in the 
quasi-particle weight, that is a possible artifact of the method 
rather than a real discontinuity of the transition. 



in favor of one of the two sublattices and the order pa- 
rameter 8 roughly measures the tendency of the system to 
form charge ordered phase. All the following results have 
been obtained for a fixed value of the local correlation 
U = 2, corresponding to a region of the phase-diagram 
with narrowing intermediate CDW-M phase separating 
the Wigner-Mott insulator from the homogeneous metal. 

Results from SB4 method are presented in Fig. 5. The 
figure shows, at V = V c i, the emergence of a charge or- 
dering instability which emerges before the Fermi liquid 
state is destroyed. At T = the slave-bosons solution 
displays the conventional critical square-root behavior 
for the order parameter, i.e. 8 = (V — Vd) 1 / 2 . At fi- 
nite temperature, however, a small "jump" in the order 
parameter can be observed, that is believed to be an arti- 
fact of the SB4 method. In the following we shall clarify 
the character of the transition using numerically exact 
methods, namely the CTQMC and DMRG at finite and 
zero temperature, respectively. 

The DMRG results for the charge ordering at zero tem- 
perature are shown in Fig. 6(a). We observe a sharp in- 
creasing of the order parameter for V ~ V c i with a typical 
squarc-root-like critical behavior for higher values of the 
inter-site correlation. This result is in good agreement 
with that obtained with SB4 method at T = 0, although 
the finite size of the DMRG effective problem does not 
permit to exactly locate the critical value of V. On the 
other hand, no trace of any discontinuoities for the order 
parameter has been observed at finite temperature, us- 
ing CTQMC (cf. Fig. 6(b)) and finite temperature Exact 
Diagonalization method (not shown). Thus, our results 
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Figure 8: (Color online) Quasi-particle residues Za,b behav- 
ior as a function of V and U = 2 at (a) T = 0.01 and (b) zero 
temperature. The figures show the linear vanishing of the 
order parameter of the Wigner-Mott transition, independent 
of the value of the temperature and expressing the continous 
character of the transition. 



put strong evidences for the continuous character of the 
charge-ordering transition, irrespective of the value of the 
temperature. 



D. Critical behavior 

Finally we have investigated the destruction of the 
charge-ordered metallic state in favor of the Wigner-Mott 
insulator, driven by the increasing inter-site correlation 
V and for large values of the local interaction U. As V 
increases, the effective hopping amplitude at the non- 
empty sublattice t e ff = t 2 /V decreases. Thus, it is 
reasonable to expect that for U ~ t e ff a Wigner-Mott 
transition takes place. Our results, obtained with differ- 
ent methods, substantiate this qualitative picture. Wc 
observe in fact a Wigner-Mott transition for V C 2(U) ~ 







Figure 9: (Color online) Order parameter 8 — ua — ub as 
a function of the non-local interaction V . Data are obtain 
from CTQMC calculations at T = 0.01 and (7 = 3. The fig- 
ure illustrates the first order-character of the direct transition 
between the homogeneous metal and the Wigner-Mott insu- 
lator, presenting evidences for the coexistence of the metallic 
and insulating phases in a narrow interval around V C 2- 



t*ff/U, that is compatible with the observed 1/U be- 
havior for the CDW metal to insulator transition line 
in the phase-diagram (cf. Sec. II). In the following we 
present a characterization of the Wigner-Mott transition 
in the EHM in terms of the vanishing of the renormal- 
izaton constants Za,b- 

Results obtained with SB4 method are presented in 
figure Fig. 7. At T = the quasiparticlc weight for the 
nearly half-filled sublattice, say A, is found to have a 
linearly vanishing behavior Za — (V C 2 — V), similar to 
that observed in the conventional Mott scenario. At fi- 
nite temperature the SB4 shows again the breakdown of 
the solution, leading to unphysical jump in the renor- 
malization constant behavior. This issue has been clar- 
ified using the more sophisticated CTQMC method at 
finite temperature. Results from CTQMC calculations 
at low temperature are shown in figure Fig. 8(a). The 
renormalization constant for the non-empty sublattice is 
found to have a linear vanishing behavior as observed 
within the SB4 solution at T = 0, thus indicating the 
continuous character of the Wigner-Mott transition. A 
very similar behavior is observed in the zero temperature 
model calculations performed with DMRG and presented 
in Fig. 8(b). 

The continuous transition from the correlated CDW- 
Mctallic phase with m* > 1 to the Mott-Wigner insula- 
tor, observed at temperatures T < T c ~ 0.01 and local 
interaction large enough U > 2, is replaced at higher tem- 
perature T > T c by a sharp first-order transition from the 
homogeneous metal with m* ~ 1 directly to the CDW- 
Insulating state. This is compatible with the shrinking 
and disappearance of the intermediate correlated CDW- 
Mctallic state with increasing temperature. This effect 



10 



is well illustrated in Fig. 9, presenting the tiny hysteresis 
cycle of the order parameter 6 — ua — tib for a large value 
of the correlation U = 3 and T = 0.01. The two curves in 
the figures are obtained following the metallic and the in- 
sulating solutions, respectively. As we can clearly see, the 
two solutions have a small region of coexistence around 
the transition point at V — V C 2 , indicating the first-order 
character of the Wigncr-Mott transition at T > T c . 

Interestingly, no trace of first-order transition have 
been found in the less correlated regime U < 2, which 
is more relevant for the interpretation of the experimen- 
tal results. In this regime the correlated quasiparticles 
arc gradually destroyed by the thermal fluctuations at 
T > T* ~ Z, leading to a thermal metal to insulator 
crossover with an associated resistivity maximum, but 
without phase separation or any first-order transition, 
compatible with out interpretation in terms of general- 
ized Pomeranchuk effect. 



IV. CONCLUSIONS AND PERSPECTIVES 

In this work we investigated the quarter filled extended 
Hubbard model, in order to describe the essential features 
of the Wigner-Mott metal-insulator transitions. Using 
single-site DMFT theory, we obtained the U-V phase- 
diagram, and carefully studied its evolution as a func- 
tion of temperature. For this model, a metallic charge 
ordered phase (CDW-M) is generally found at low tem- 
peratures, separating a Wigner-Mott insulator at strong 
coupling from a homogeneous Fermi liquid at weak cou- 
pling. This intermediate CDW-M phase is the one show- 
ing the most interesting features, chiefly the emergence 
of strong correlation effects signaled by large effective 
mass m* enhancements. Physically, this reflects the pres- 
ence of heavy quasiparticles existing only below a char- 
acteristic energy scale E* ~ Tp/m*, which vanishes at 
the metal-insulator transition. Indeed, we predict that 
such a correlated metallic state can easily be suppressed 
either by increasing the temperature beyond a modest 
temperature T* ~ E* , or by applying modest polariza- 
tion fields H* ~ T* - in striking agreement with the 
experiments . In addition, we demonstrated that the re- 
gion occupied by the CDW-M phase shrinks as a function 
of increasing temperature, thus producing an interesting 
Pomeranchuk- like effect. We presented an interpretation 



of this effect in terms of the entropic destruction of the 
strongly renormalized Fermi liquid (CDW-M) in favor 
the Wigncr-Mott insulator, similar to the early idea pro- 
posed by Pomeranchuk in the context of the 3 He solidi- 
fication. 

The extended Hubbard model we used certainly cannot 
be regarded as a realistic or quantitatively accurate rep- 
resetation of the 2DEG materials. Ours is an approach 
complementary to that provided by first principle (e.g. 
diffusion Monte Carlo) studies of realistic models 49 of 
2DEG - yet with surprisingly similar results. Both ap- 
proaches portray a picture of a strongly correlated elec- 
tron fluid featuring a single characteristic energy scale. 
But precisely by virtue of its simplicity, our model calcu- 
lation makes it possible to unravel the mechanism behind 
the puzzling behavior in the metal-insulator transition re- 
gion. It demonstrates how the tendency to charge order- 
ing reinforces the transmutation of conduction electrons 
into local magnetic moments - a fundamental physical 
process behind the phenomenon Wigner-Mott localiza- 
tion. 

Many quantitative aspects of our model can and should 
be improved. In particular, our lattice model does not 
do justice to dynamical charge fluctuations, which should 
further enhance the role of Coulomb correlations even 
in absence of long-range charge order. Indeed, results 
supporting the robustness of the intermediate correlated 
metallic phase, have been recently obtained in Ref. [50] 
using complementary methods which emphasized the im- 
portant role of longer-range Coulomb interactions. In 
addition, the influence of disorder also needs to be ad- 
dressed in the context of the Wigner-Mott scenario we 
propose. These effects can be naturally incorporated in 
our framework using the recently developed "extended" 
DMFT approaches 01,52 . This interesting and important 
task opens an interesting avenue for future work. 
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